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Abstract 

Background 

A highly pathogenic human coronavirus (CoV), Middle East respiratory syndrome coronavirus 
(MERS-CoV), has emerged in Jeddah and other places in Saudi Arabia, and has quickly 
spread to European and Asian countries since September 2012. Up to the 1 st October 2015 it 
has infected at least 1593 people with a global fatality rate of about 35%. Studies to understand 
the virus are necessary and urgent. In the present study, MERS-CoV main protease (M pro ) is 
expressed; the dimerization of the protein and its relationship to catalysis are investigated. 

Methods and Results 

The crystal structure of MERS-CoV M pro indicates that it shares a similar scaffold to that of 
other coronaviral M pro and consists of chymotrypsin-like domains I and II and a helical 
domain III of five helices. Analytical ultracentrifugation analysis demonstrated that MERS- 
CoV M pro undergoes a monomer to dimer conversion in the presence of a peptide substrate. 
Glul 69 is a key residue and plays a dual role in both dimerization and catalysis. The muta¬ 
genesis of other residues found on the dimerization interface indicate that dimerization of 
MERS-CoV M pro is required for its catalytic activity. One mutation, M298R, resulted in a sta¬ 
ble dimer with a higher level of proteolytic activity than the wild-type enzyme. 

Conclusions 

MERS-CoV M pro shows substrate-induced dimerization and potent proteolytic activity. A 
critical assessment of the residues important to these processes provides insights into the 
correlation between dimerization and catalysis within the coronaviral M pro family. 


Introduction 

A highly pathogenic human coronavirus (CoV) 1 , Middle East respiratory syndrome coronavi¬ 
rus (MERS-CoV), emerged in Jeddah and other places in Saudi Arabia in September 2012 and 
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spread to some European, African and Asian countries in recent years [1-3]. The virus causes 
symptoms similar to severe acute respiratory syndrome coronavirus (SARS-CoV), but an addi¬ 
tional component is also involved, namely acute renal failure [4], Up to the 1 st October 2015, 
1593 people had been infected with MERS and this has led to 568 reported deaths globally 
(World Health Organization, global alert and response, http://www.who.int/csr/don/ 
01-october-2015-mers-jordan/en/). A serological study of major livestock in Saudi Arabia sug¬ 
gested dromedary camels may be an original reservoir [5, 6]; although recent studies have 
identified bats may also as the suspected host as a number of bat Co Vs show high sequence 
similarity to SARS-CoV and MERS-CoV [7, 8] and a bat-CoV was discovered that readily 
infect human cells using ACE2 as the receptor [9], Nevertheless, human-to-human transmis¬ 
sion of MERS-CoV has been confirmed. This May, an infected traveler from the Middle East 
region returned to his home country of the Republic of Korea and caused an outbreak that 
revolved around health facilities (http://www.who.int/csr/don/07-july-2015-mers-korea/en/). 
These findings indicate that it is possible for the virus to spread globally and as such it poses a 
significant threat to world health and the world economy in general. Therefore studies that 
help to understand this virus and aid the development of antiviral drugs or other therapies are 
important. 

Similar to other Co Vs, after entering the host cells, the nonstructural polyproteins (ppla 
and pplab) of MERS-CoV are synthesized and then cleaved by two coronaviral proteases, a 
main protease (M pro ) (EC 3.4.22.69) and a papain-like protease (EC 3.4.22.46) [10]. This cleav¬ 
age is considered to be a leading process which is required for viral maturation [11-14]. The 
MERS-CoV M pro , namely nsp5 of the ppla proteins (residue 3248-3553), has been identified 
[15]. Like other M pro , there is a catalytic dyad that consists of a His residue and a Cys residue 
[16-20]. Sequence alignment suggests that MERS-CoV M pro , in a similar manner to other 
known M pro , has a chymotrypsin-like architecture consisting of a catalytic core (domain I and 
II) and a helical domain III; sequence identities of CoV M pro protein range from 50% to 80% 
(SI Fig). Recently the crystal structure of an inactive MERS-CoV M pro C148A mutant has con¬ 
firmed this similarity and the results also suggest that this protein forms a dimer [15]. Further¬ 
more, based on the identification of eleven canonical cleavage sites, the MERS-CoV M pro 
should be able to recognize and cleave at the (L/M)-Q-j-(A/S) conserved sequence, which is 
essential for most CoV M pro -mediated processing [21, 22]. However, up to now, the correlation 
between the protein’s structure and the catalytic process remains unclear. 

In the present study we expressed and purified the MERS-CoV M pro using the authentic N- 
terminus via an Escherichia coli system. The crystal structure of the MERS-CoV M pro at 3.0-A 
resolution is reported. The quaternary structural changes in the MERS-CoV M pro in the 
absence and presence of peptide substrates were investigated by analytical ultracentrifugation 
(AUC). The results of kinetic activity assays indicated that MERS-CoV M pro exhibits potent 
proteolytic activity that is associated with a pattern of cooperativity. Some critical residues for 
dimerization of the protein and catalysis by the protein were verified by site-directed mutagen¬ 
esis. The present studies provide a foundation for an understanding of the mechanism that 
controls the monomer-dimer switch at work in MERS-CoV M pro . 

Materials and Methods 

Expression Plasmid Construction 

The sequence of the MERS-CoV M pro (GenBank accession number AHC74086; polyprotein 
residues 3248-3553) was synthesized (MDBio Inc.), digested by Ndel-Xhol and then inserted 
into the vector pET-28a(+) (Novagen). In this construct, the 6 x His tag is retained at the 
N-terminus. To remove the fusion tag and generate an authentic N-terminus for protein 
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purification, the codons of the thrombin cutting recognition sequence and a Nde I cutting site 
were removed and then inserted the codons of Leu-Arg-Leu-Lys-Gly-Gly into the above vector. 
The forward primer sequence for site-directed mutagenesis was 5’-CATCACAGCAGCGGCCT 
GCGTCTGAAAGGCGGCAGCGGTTTGGTGAAAATG-3’ and the reverse primer was 5’-CATTTT 
CACCAAACCGCTGCCGCCTTTC AGACGCAGGCCGCTGCTGTGATG-3’. The reading frame of 
the final plasmid was confirmed by sequencing. 

Expression and Purification of MERS-CoV M pro 

The expression vector was transformed into E. coli BL21 (DE3) cells (Novagen). Cultures were 
grown in 0.8 liters of LB medium at 37°C for 4 h, induced with 0.4 mM isopropyl-(3- D -thioga- 
lactopyranoside, and then incubated overnight at 20°C. After centrifuging at 6,000 x g at 4°C 
for 15 min, the cell pellets were resuspended in lysis buffer (20 mM Tris, pH 8.5, 250 mM 
NaCl, 5% glycerol, 0.2% Triton X-100, and 2 mM fS-mercaptoethanol) and then lysed by soni- 
cation. The crude extract was then centrifuged at 12,000 x g at 4°C for 25 min to remove the 
insoluble pellet. Next the supernatant was incubated with 1-ml Ni-NTA beads at 4°C for 1 h 
and then loaded onto an empty column. After allowing the supernatant to flow through, the 
beads were washed with washing buffer (20 mM Tris, pH 8.5, 250 mM NaCl, 8 mM imidazole, 
and 2 mM P-mercaptoethanol). The SARS-CoV papain-like protease [12] (1 mg in 100 mM 
phosphate buffer (pH 6.5)) was then added and incubated for 3 h. The SARS-CoV papain-like 
protease digestion, which removed the 6 x His tag and Leu-Arg-Leu-Lys-Gly-Gly fragment, 
resulted in a native protein product with an authentic N-terminus. The digest was allowed to 
flow through and then loaded onto a S-100 gel-filtration column (GE Healthcare) equilibrated 
with running buffer (20 mM Tris, pH 8.5,100 mM NaCl, and 2 mM dithiothreitol). The purity 
of the fractions collected was analyzed by SDS-PAGE and the protein was concentrated to 30 
mg/ml by Amicon Ultra-4 10-kDa centrifugal filter (Millipore). 

Protein Crystallography 

Crystals of the MERS-CoV M pro were obtained at 295 K by the sitting-drop vapor-diffusion 
method. The protein solution was set up at 5 mg/ml and the reservoir solution consisted of 0.1 
M Tris, pH 8.4,15% (w/v) PEG 4000 and 0.2 M sodium acetate. Clusters of needle crystals 
appeared in 2 days and were used for micro-seeding. Single cystals of rectangle shape and with 
dimensions of 0.3-0.5 mm were obtained in less than a week. All crystals were cryoprotected in 
the reservoir solution with 15% glycerol and were flash-cooled in liquid nitrogen. 

Data collection, structure determination and refinement 

X-ray diffraction data were collected at 100 K on the SPXF beamline 13C1 at the National Syn¬ 
chrotron Radiation Research Center, Taiwan, ROC, using a ADSC Quantum-315r CCD detec¬ 
tor (X-ray wavelength of 0.976 A). The diffraction images were processed and scaled using the 
HKL-2000 package [23]. The structure was solved by the molecular replacement method by 
Phaser [24] using the structure of SARS-CoV M pro R298A mutant (PDB entry 4hi3; [25]) as 
the search model. Manual rebuilding of the structure model was performed using Coot [26]. 
Structure refinement was carried out using REFMAC [27]. The data-processing and refinement 
has been deposited in the Protein Data Bank (PDB entry 5c3n). 

Steady-State Kinetic Analysis 

The colorimetry-based peptide substrate, TSAVLQ-para-nitroanilide (TQ6-pNA) (purity 95- 
99% by HPLC; GL Biochem Ltd, Shanghai, China), was used to measure the proteolytic activity 
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of MERS-CoV M pro and its mutants throughout the course of the study as described previously 
[25, 28] . This substrate is cleaved at the Gln-pNA bond to release free pNA, resulting in an 
increase in absorbance at 405 nm. The absorbance at 405 nm was continuously monitored 
using a Jasco V-550 UV/VIS spectrophotometer. The protease activity assay was performed in 
10 mM phosphate (pH 7.6) at 30°C. The substrate stock solution was 1600 pM and the working 
concentrations were from 25 to 1200 pM. In the substrate titration assay, the concentration of 
MERS-CoV M pro and its mutants, V4R, T126S, E169A, M298R and T126S/M298R was 0.3, 

0.4, 0.7,1.2, 0.15 and 0.26 pM, respectively, while that of SARS-CoV M pr ° was 1.1 pM. Steady 
state enzyme kinetic parameters were obtained by fitting the initial velocity (v 0 ) data to the 
Michaelis-Menten Eq (1) 


_ UE][S] 
K m + [S] 


( 1 ) 


where k^ is the catalytic constant, [E] is the enzyme concentration, [S] is the substrate concen¬ 
tration and K m is the Michaelis constant of the substrate. The program SigmaPlot (Systat Soft¬ 
ware, Inc., Richmond, CA) was used for the data analysis. 

To assess the cooperativity effect, the kinetic parameters were obtained by fitting the initial 
velocities to the Hill Eq (2) 


_ UE][Sf 

K' + [Sf 


( 2 ) 


where K’ is a constant that is related to the dissociation constant and h is the Hill constant. 


Analytical ultracentrifugation analysis 

AUC was performed on a XL-A analytical ultracentrifuge (Beckman Coulter) using an An-50 
Ti rotor [11,12, 25, 28-30]. The sedimentation velocity experiments were carried out using a 
double-sector epon charcoal-filled centerpiece at 20°C with a rotor speed of 42,000 rpm. Pro¬ 
tein solutions of 0.05 to 0.5 mg/ml (330 pi) and reference (370 pi) solutions, both containing 
D 2 0, were loaded into the centerpiece. The absorbance at 280 nm was monitored in a continu¬ 
ous mode with a time interval of 300 s and a step size of 0.003 cm. Multiple scans at different 
time intervals were then fitted to a continuous c(s) distribution model using the SEDFIT pro¬ 
gram [31]. Additionally, the results with the various different protein concentrations were 
globally fitted to a monomer-dimer self-association model using the SEDPHAT program to 
calculate the dissociation constant (Kj) [32]. 

To measure the substrate-induced dimerization, the active enzyme centrifugation (AEC) 
[33] was performed. Briefly, MERS-CoV M pro of 15 pi (1 mg/ml) was added into the small well 
of the band-forming centerpiece before the cell assembled. Then 330 pi of peptide substrate at 
0, 200 and 400 pM in D 2 0 were respectively loaded into the bulk sample sector space. At a 
rotor speed of 42,000 rpm, the protein solution flowed into the substrate-containing channel 
and form a protein band. It can be detected by absorbance at 250 nm. During the centrifuga¬ 
tion, the sediment protein continuously met and cleaved the substrate, which can be detected 
by absorbance change at 405 nm. The dataset from the multiple scans at 250 nm at various 
time intervals were fitted to a continuous c(s) distribution model using the SEDFIT program 
[31], while the first five scans (0-30 min) at 405 nm were used to derive the product concentra¬ 
tion and then initial velocity values. 
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Isothermal titration calorimetry (ITC) 

The protocol followed that of previous studies [28] with some modifications. Apparent dissoci¬ 
ation constants and stoichiometry of the enzyme-ligand interactions were measured by a Ther¬ 
mal Activity Monitor 2277 from TA instruments (New Castle, DE). Calorimetric titrations of 
the peptide substrate TQ6-pNA (0.5 mM in a 250-pi syringe) and M pro (6 pM in a 4-ml 
ampoule) were carried out at 25°C in 10 mM phosphate buffer (pH 7.6). The peptides were 
titrated into the enzyme using a 10-pl aliquot for each injection with a time interval of 20 min. 
A control experiment in the absence of enzyme was performed in parallel to correct for the 
dilution of heat. The data obtained was then analyzed by integrating the heat effects normalized 
against the amount of injected protein using curve-fitting based on a 1:1 binding model. This 
involved the use of Digitam software (TA instruments, New Castle, DE). 


Results and Discussion 

Recombinant MERS-CoV M pro preparation 

As part of the present study, an expression vector was constructed and the BL21 (DE3) STAR 
(Invitrogen) strain of E. coli were used to express MERS-CoV M pro . Unlike SARS-CoV M pro 
[25, 28], the MERS-CoV M pro with 6 x His-tag retained at the C-terminus cannot be expressed. 
Instead, the bacteria are able to express the M pro when there is a N-terminal 6 x His-tag fusion 
that can be removed during the purification. However, thrombin digestion leaves two extra res¬ 
idues (Gly-Ser) at the N-terminus of M pro , resulting in protein with no proteolytic activity 
(data not shown). Therefore we used SARS-CoV papain-like protease [12, 30, 34], which is a 
highly active viral deubiquitinase and does not leave any residues at the N-terminus of M pro . 
After gel-filtration, the purity of authentic N-terminus M pro was about 99% (S2 Fig). The size 
of the MERS-CoV M pro was found to be close to 30 kDa, while any uncut protein was located 
at higher molecular weight position. The typical yield was about 10 mg after purification from 
0.8 liter of E. coli culture. 


Overall structure of MERS-CoV M pro 

The structure of the MERS-CoV M pro was determined at 3.0 A resolution by X-ray crystallog¬ 
raphy (Table 1 and Fig 1 A). The crystal packing belonged to space group C222 1; with unit-cell 
parameters a = 87.2, b = 94.0, c = 155.1 A and a = [1 = y = 90°. The final atomic model contain¬ 
ing two M pro molecules in a crystallographic asymmetric unit agrees well with the crystallo¬ 
graphic data and the expected values of geometric parameters (Table 1). There are no residues 
in the disallowed region of the Ramachandran plot, while 81.3% of the residues are in the most 
favored region. 

The overall dimeric structure of M pro is similar to that of SARS-CoV M pr °; although a rela¬ 
tive shift of 10° to 30° could be observed for the two domain III within the dimer (S3A Fig). 
Indeed, the r.m.s. distance between equivalent Ca atoms of the domain I+II of the two struc¬ 
tures is 0.9 A, while that between the domain III of the two structures is 3.1 A. Compared with 
the structures of the ligand-bound complex (PDB entry 4YLU) and C148A mutant (PDB entry 
4WME) [14, 15], the r. m. s. distance is 0.8 and 0.7 A over 540 Ca atom pairs, respectively (S3B 
Fig). This indicates that the dimeric structures show no significant difference; although the 
present structure is a free enzyme and the other two structures involve enzyme-ligand com¬ 
plexes and higher resolution. Besides, there is minor difference between the present structure 
and bat-CoV HKU4 M pr ° (PDB entry 2YNA), with 80% sequence identity, as the r. m. s. dis¬ 
tance over 539 Ca atom pairs is 0.8 A (S3B Fig). 
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Table 1. Summary of crystallographic information for MERS-CoV M pro . 
Data Collection 


Space group 

C222t 

Cell dimensions 

a, b, c (A) 

87.2, 94.0, 155.1 

a, 3, Y O 

90, 90, 90 

Resolution 3 (A) 

30-3.0 (3.11-3.0) 

Ft merge b (%) 

17 (62.9) 

/ / ol 

10.2 (3.1) 

Completeness (%) 

94.6 (95.2) 

Redundancy 

5.5 (5.6) 

Refinement 

Number of reflections 

11,675 (1,549) 

R factor^ (%) 

21.8 

Free R factor d (%) 

28.2 

Number of atoms 

Protein 

4,580 

Water 

0 

Average B-factors for protein atoms (A 2 ) 

81.6 

R.m.s. deviations 

Bond length (A) 

0.009 

Bond angles (°) 

1.5 

Ramachandran plot statistics (%) 

Most favored region 

81.3 

Additional allowed region 

13.6 

Generously allowed region 

5.2 

Disallowed region 

0 


a The numbers in parentheses are for the highest-resolution shell. 

b fi mom = 22^2 \l h i~ ( l h ) |/^2 22hi’ where hi is the integrated intensity of a given reflection and (l h ) is the 

hi hi 

mean intensity of multiple corresponding symmetry-related reflections. 

C R = 22 | F a h - F c h | /22 where F° and F c h are the observed and calculated structure factors, respectively. 

h h 

d Free R is R calculated using a random 5% of data excluded from the refinement. 
doi:10.1371/journal.pone.0144865.t001 

Interestingly, the dimerization interface situation with the MERS-CoV M pro was found to be 
different to that of SARS-CoV M pro where there are four amino acid pairs with intermolecular 
polar interactions (Seri.. .Glul66, Arg4.. .Glu290, Serl23.. .Arg298 and Serl39.. .Gln299). 
There are only two pairs of intermolecular hydrogen bonds, Seri.. .Glul69 and Serl42.. .Gln299 
that are associated with the dimer surface of MERS-CoV M pro according to the current structure 
(Fig IB). This led us to compare the dimerization and catalytic activity of the two types of M pro . 
In the present study, in addition to using the wild-type MERS-CoV M pro , we also mutated several 
residues at the dimerization interface in order to evaluate their role in dimerization and catalysis 
of MERS-CoV M pro (see below). 

Correlation between dimerization and catalysis of MERS-CoV M pro 

To compare catalysis between the two M pro , TQ6-pNA, a peptide substrate for SARS-CoV 
M pro [25,28], was used to measure the proteolytic activity. At first, the dependence of the initial 
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A 



Fig 1. The structure of MERS-CoV M pro . (A) The overall structure of the dimeric M pro . The two protomers 
are shown as a ribbon and as a surface model, respectively. The negatively and positively charged regions 
on the molecular surfaces are colored red and blue. The spheres indicate the catalytic dyad, His41 -Cysl 48. 
(B) Stereo view of an overlay of the dimerization interface and the active site of MERS-CoV M pro (in cyan and 
orange) with that of SARS-CoV M pro (grey; [16]). The red dashed lines indicate polar interactions between the 
two protomers of MERS-CoV M pro , while the black dashed lines show polar interactions between those of 
SARS-CoV M pro . The C atoms of the modeled substrate P4-P1 residues (from the structure of Cl 48A mutant 
[15]) are colored magenta. The structural figures in this paper were produced using PyMol (http://www.pymo!. 

org/). 

doi:10.1371/journal, pone. 0144865.g001 

velocity on enzyme concentration was analyzed and showed a nonlinear upward correlation 
(Fig 2A). The pattern is similar to that of SARS-CoV M pro , as the monomeric M pro may not 
have catalytic activity [28]. However, MERS-CoV M pro displayed a sigmoid curve for its rate 
constant pattern at various substrate concentrations (Fig 2B, open circles); this contrast with 
SARS-CoV M pro , which exhibited a classical saturation curve (S4 Fig). The results were then 
fitted to the Hill equation (Eq 2) in order to evaluate the kinetic parameters (Table 2). The k cat 
(2.33 s" 1 ) of MERS-CoV M pro is close to that of SARS-CoV M pro (2.11 s" 1 ), while the Hill con¬ 
stant was 1.8, suggesting a significant degree of positive cooperativity among the M pro proto¬ 
mers. The comparable activity levels of the two M pro in the present study is dissimilar to the 
results obtained during a recent study in which the activity of the MERS-CoV M pro was found 
to be 5-fold lower than that of SARS-CoV [14]. Using different substrates may cause the differ¬ 
ence. Tomar et al. [14] used a longer peptide substrate with residues present at both P and P’ 
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WT (uM) 



TQ6-pNA (uM) 


Fig 2. Proteolytic activity assay of MERS-CoV M pro and its mutants. (A) Plot of initial velocities (v,) versus 
the concentration of wild-type MERS-CoV M pro . The concentration of substrate (TQ6-pNA) was 600 pM. The 
line represented the best-fit results according to the nonlinear dependence equation [28], (B) and (C) Plots of 
rate constant (k obs ) versus the concentration of peptide substrate for MERS-CoV M pro (by circles), the V4R 
mutant (by lower triangles), theT126S mutant (by squares), the E169A (by diamonds), the M298R (by upper 
triangles) and the T126S/M298R (by hexagons) are indicated. The lines represented the best-fit results 
according to the Hill equation (Eq 2; wild-type, V4R, T126S and E169A mutants) or the Michaelis-Menten 
equation (Eq 1; M298R and T126S/M298R mutants). The protein concentrations of the wild-type, V4R, 
T126S, E169A, M298R and T126S/M298R mutants used for the assay were 0.3, 0.4, 0.7,1.2, 0.15 and 
0.26 pM, respectively. All the assays were performed in 10 mM phosphate (pH7.6) at 30°C and repeated 
twice to ensure reproducibility and the error bars were shown. The kinetic parameters are shown in Table 2. 

doi:10.1371/journal, pone.0144865.g002 
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PLOS 


Table 2. The kinetic parameters and dissociation constants of MERS-CoV M pro . 


Proteins 


Kinetic Parameters 3 


Dissociation Constant* 1 

K m (uM) 

kcat (s' 1 ) 

h 

No substrate(pM) 

With 600 pM substrate (pM) 

MERS-CoV M pro 

Wild-type 

- 

2.33 ±0.13 

1.8 ± 0.04 

7.7 ±0.3 

0.7 ±0.04 

V4R 

- 

0.96 ± 0.05 

2.7 ± 0.2 

23.0 ± 0.4 

15.2 ±0.3 

T126S 

- 

0.56 ± 0.04 

2.0 ± 0.2 

33.7 ±0.9 

13.9 ±0.1 

E169A 

- 

0.41 ± 0.02 

2.1 ± 0.1 

14.3 ±0.2 

14.1 ± 0.5 

M298R 

181.0 ±24.0 

7.91 ± 0.49 

- 

1.1 ±0.1 

0.7 ±0.01 

T126S/M298R 

419.4 ±63.9 

4.63 ± 0.37 

- 

2.8 ±0.1 

0.9 ±0.01 

SARS-CoV Mpro 

890 ±130 

2.11 ± 0.15 

- 

0.7 ± 0.02 

1.7 ± 0.03 c 


3 Kinetic data of SARS-CoV M pro and MERS-CoV M298R and T126S/M298R mutants were fitted to the Michaelis-Menten equation (Eq 1), while those of 
the others were fitted to the Hill equation (Eq 2). The R sqr were from 0.985 to 0.999, respectively. All the assays were repeated twice and the average 
values were used for the fitting. 

b The values were derived from a global fit of the AUC data to a monomer-dimer self-association model by SEDPHAT [32]. The experiments for the assay 
were obtained at protein concentration of 1.5 to 30 pM. 
c The value was from our previous studies for comparison [28]. 


doi: 10.1371 /jou rnal. pone.0144865 ,t002 


site. However, in the present studies we used a peptide substrate that contains only P site resi¬ 
dues. Besides, they utilized FRET substrate and could only be used at low substrate concentra¬ 
tions to prevent the inner-filter effect, while we are able to use higher substrate concentrations 
to capture the kinetic parameters. 

The cooperativity phenomenon associated with the MERS-CoV M pro is similar to that of 
the SARS-CoV M pro R298A/L monomer mutants; these were found to show monomer to 
dimer conversion during catalysis [28] . As a result of the above, we investigated the quaternary 
structure of the M pro by AUC (Fig 3). The cumulative spectra (Fig 3A) were analyzed using the 
continuous c(s) distribution model and the results suggested that MERS-CoV M pro is a mono¬ 
mer in phosphate buffer (S5 Fig) and this contrasts with a distribution of 30% monomer and 
70% dimer in the presence of 600 pM TQ6-pNA (Fig 3B). We also measured the size distribu¬ 
tion of M pro at various TQ6-pNA concentrations (S6 Fig). The results indicated that the sedi¬ 
mentation coefficient of the major species was shifted as the substrate dosage changed (S6B 
Fig). More substrate led to the major species moving close to the dimer position. 

However, before the centrifugation, the enzyme had been mixed with the substrate and the 
catalysis began, resulting in a mixture of substrate and product with enzyme. It is unable to 
confirm that our observation is a substrate-induced or substrate/product-induced dimeriza¬ 
tion. To solve this, a modified AUC technique, AEC [33], was utilized to detect the quaternary 
structure change in the absence and presence of substrate (Fig 4). Here the enzyme solution 
was put into the small well of band-forming centerpiece and then flowed into the substrate- 
containing channel when the centrifugation began. During the centrifugation, the protein layer 
gradually sediment and continuously met peptide substrates. Not surprisingly, there was a 
broad distribution between the monomer and dimer species in the presence of 200 pM peptide 
substrate, while a major species shifted to the dimer in 400 pM substrate (Fig 4B). It suggests 
that MERS-CoV M pro acts as a rapid self-associated and substrate-induced dimerization. Using 
different strategies, Tomar et al. [14] confirmed that inhibitor binding can also induce and 
maintain the dimerization of M pro . On the other hand, we measured the velocity of the product 
formation during the centrifugation; although the rate (0.017 pM/s) is 10-fold lower than that 
by the spectrometric assay (Fig 4B). 
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Fig 3. The continuous size distribution change of MERS-CoV M pro . (A) Typical trace of absorbance at 
250 nm of the enzyme during the SV experiment. The protein concentration was 0.5 mg/ml. For clarity, only 
every fifth scan is shown. The symbols represent the experimental data and the lines represent the results 
fitted to the Lamm equation using the SEDFIT program [31]. (B-G) Continuous c(s) distribution of wild-type, 
V4R, T126S, El 69A, M298R and T126S/M298R mutants. The distributions in D 2 0 containing 10 mM 
phosphate (pH 7.6) are shown by solid lines and those in the same buffer but with 600 pM TQ6-pNA 
substrate are shown by dashed lines. The left vertical dotted line indicates the monomer position and the right 
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dotted line represents the dimer position. The residual bitmap of the raw data and the best-fit results are 
shown in the insets. 

doi:10.1371/journal, pone. 0144865.g003 

To quantitatively characterize the monomer-dimer equilibrium of M pro in the absence and 
presence of substrate, the AUC results at protein concentrations of 1.5, 6 and 15 pM were glob¬ 
ally fitted to the monomer-dimer self-association model (Table 2). The analysis indicated that 
the K d value of MERS-CoV M pro in the absence of substrate was 7.7 pM and this decreased 
11-fold in the presence of substrate, which brings it close to the value for SARS-CoV M pro 
(1.7 pM; Table 2). Thus it can be concluded that, like the SARS-CoV M pro R298A mutant [25], 
the presence of substrate is able to induce the dimerization of MERS-CoV M pro . In addition, 
although the K d values of wild-type SARS-CoV M pro without or with substrates show no signif¬ 
icant difference (Table 2), it was possible to detect substrate-induced dimerization at a protein 
concentration of 1 pM by AEC [33]. 

Previous studies have demonstrated that a conserved residue, Glul66, plays a pivotal role 
in connecting the substrate binding site with the dimerization interface of SARS-CoV M pro 
[28]. Here the equivalent residue, Glul69, was mutated and its influence on MERS-CoV M pro 
evaluated (Fig IB). Not surprisingly, compared to the wild-type enzyme, the activity (k cat ) of 
E169A showed a 6-fold decrease (Fig 2B, open diamonds and Table 2). Furthermore, AUC 
analysis suggested that this enzyme consisted of a single species close to monomeric form 
even in the presence of substrate (Fig 3E). The K d values of the E169A mutant without sub¬ 
strate was 14 pM, which is 2-fold higher than that of the wild-type and this did not decrease in 
the presence of substrate (Table 2). The results suggest that mutation of Glul69 is able to 
block substrate-induced dimerization and that this results in a decrease in enzyme activity. 
Unexpectedly, as well as the monomer, some octamer (6.3%) with a sedimentation coefficient 
of 4.9 S was observed in the presence of substrate (Fig 3E). Previous studies have suggested 
that a super-active octamer of SARS-CoV M pro can be locked by 3D domain swapping [35]. 

In this study, the presence of the octamer form of the MERS-CoV M pro El69A mutant in the 
presence of substrate may explain why this mutant is not totally inactive. Taking the above as 
a whole, the dual role of the conserved Glu residue in catalysis and dimerization is consistent 
for both M pro . 

Our results confirmed that there are fewer intermolecular polar interactions at the dimeriza¬ 
tion interface of MERS-CoV M pro and this results in the enzyme being in the monomer form 
in aqueous buffer; this contrasts with SARS-CoV M pro , which is mostly in the dimer form in 
similar circumstances due to the greater number of intermolecular interactions (Figs IB and 
3B). Based on the sequence alignment, three residues in the dimerization interface vary in coro- 
naviral M pro sequences (SI Fig). In SARS-CoV M pro , residues Arg4, Serl23 and Arg298 are 
different to the equivalent ones in MERS-CoV M pro , which are Val4, Thrl26 and Met298, 
respectively. Based on the above, three single mutants, V4R, T126S and M298R, were generated 
and their proteolytic activity and dimerization were assessed. Unexpectedly, both the V4R and 
the T126S mutant showed a higher K d for the dimer to monomer and a lower level of activity 
than wild-type M pro (Fig 2B and Table 2). K d values were decreased by 1.5-fold to 2.4-fold for 
the two mutants in the presence of substrate, which suggests reduced substrate-induced dimer¬ 
ization (Fig 3C and 3D). Based on the current structure, residue Val4 showed a hydrophobic 
contact with another protomer’s residue Glyl41 (Fig IB). Mutation of valine to arginine may 
lose the contact. Furthermore, the side chain of residue Vall31 of MERS-CoV M pro is hydro- 
phobic while that of the equivalent residue, Cysl28 of SARS-CoV M pro , is hydrophilic (Fig IB). 
Vall31 is close to the Arg4 and will disfavor the electrostatic interaction of Arg4.. .Glu290 in 
MERS-CoV M pro . These variance may result in V4R failed to form a stable dimer. For the 


PLOSONE | DOI:10.1371/journal.pone.0144865 December 14, 2015 


11/18 





PLOS 


ONE 


Dimerization and Catalysis of MERS-CoV Main Protease 



< 6.3 6.4 6.5 6.6 6.7 6.8 6.9 

Radius (cm) 



Fig 4. Charactering the substrate-induced dimerization of MERS-CoV M pro by AEC [33]. (A) A trace of 
absorbance at 250 nm for the enzyme in the presence of 200 pM of TQ6-pNA during the experiment. For 
clarity, only every twice scan and every fifth measuring points per scan are shown. The protein of 15 pg in 10 
mM phosphate (pH7.6) was used. The substrate was dissolved in D 2 0 to give a higher density, which can 
maintain the protein band when the centrifugation begins. The symbols are experimental data and the lines 
are the results fitted to the Lamn equation using SEDFIT program [31]. The continuous size distribution of the 
best-fit result was shown in panel B colored by green. (B) Continuous c(s) distributions of MERS-CoV M pro at 
the concentration of TQ6-pNA of 0, 200 and 400 pM. The labels M and D showed the position of the monomer 
and dimer species. (C) Monitoring the enzyme activity of MERS-CoV M pro during the centrifugation. The 
same cell in panel A was followed and the absorbance at 405 nm trace for the released product (pNA) was 
shown. We used the spectra of the first 30 min to calculate the product concentration at different time. The 
inset plot showed the product versus time and the slope of the line represented the initial velocity. 

doi:10.1371/journal, pone. 0144865.g004 
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T126S mutant, after compared with the other two MERS-CoV M pro structures (PDB entry 
4YLU and 4WME), we found that the side chain of Thrl26 is free of rotation. It is able to make 
a hydrophobic interaction with the residue Tyrl21, resulting in a 180°-rotation of the phenol 
ring, and lead to a hydrogen-bond with the backbone amide of Leul44 (Fig IB). It further 
results in the side chain of Leul44 toward another protomer’s Ile300 and make a hydrophobic 
contact for the two protomers. Mutation of Thrl26 to serine will lose this contact and may dis¬ 
favor the dimerization. 

By way of contrast, the M298R mutant resulted in a stable dimer form in phosphate buffer 
(Fig 3F). The K d value for the mutant in the absence and presence of substrate were 1.1 and 
0.7 pM, respectively, which are very close to those for SARS-CoV M pro (Table 2). Furthermore, 
the stable dimer form showed higher proteolytic activity and the mutation transformed the 
enzymes rate constant pattern at various substrate concentrations into a classical saturation 
curve (Fig 2C, open triangles). In addition, the K m of the M298R mutant is 5.8-fold lower than 
that of SARS-CoV M pro , which suggests a higher substrate binding affinity (Table 2). It can be 
concluded that mutation of residue Met298 to arginine within the MERS-CoV M pro results in 
the stabilization of the dimer formation, which in turn gives rise to more efficient catalysis. 
Based on our structure, Arg298 is able to make a hydrogen bonding interaction with Thrl26. 
Residue Thrl26 can be replaced by serine because the T126S/M298R double mutant also 
shows similar dimerization characteristics (Fig 3G) and saturated catalytic pattern to that of 
the M298R mutant (Fig 2C, open hexagons and Table 2). However, it can only be achieved in 
the presence of Arg298, not Met298. 

Substrate binding analysis of MERS-CoV M pro by ITC 

The sigmoid nature of the curve describing the rate constant pattern at various substrate con¬ 
centrations (Fig 2B) means that it is not possible to obtain a K m value for this enzyme, which 
would allow us to evaluate the substrate binding affinity of the wild-type MERS-CoV M pro and 
its mutants; the exceptions being the M298R single mutant and the T126S/M298R double 
mutant (Fig 2). To further delineate the binding of substrate to the enzyme, ITC was used to 
measure the K d for the substrate (or substrate/product)-enzyme complex and the binding stoi¬ 
chiometry (N) (Fig 5). During the titration, the enzymatic hydrolysis might produce additional 
heat, resulting in higher AH. So we can only compared the N and K d of the wild-type M pro with 
those of E169A and M298R mutants. The three enzymes exhibited similar N (0.89 to 1.06) and 
K d (14.2 to 20.3 pM). This suggested that the monomeric and dimeric M pro show quite the 
same substrate binding affinity. Such phenomenon is also found in monomeric and dimeric 
SARS-CoV M pro , whose N and K d for the same substrate were 0.97 to 1.05 and 29.9 to 33.8 pM, 
respectively [28]. With the K m ((^+1^)/^), Iq.^ and K d (k^/ki), we are able to calculate the 
k, and k_,, the rate of the association and dissociation of the enzyme-substrate complex. The 
kj and k. : for the M298R mutant and substrate is 0.049 s _1 pM 1 and 0.98 s" 1 , while those for 
SARS-CoV M pr ° is 0.00246 s 1 pM 1 and 0.084 s' 1 . The rate constants for M298R mutant 
showed 20- and 12-fold higher than those for SARS-CoV M pro . The more rapid association 
and relatively slower dissociation of enzyme-substrate complex may be used to explain why the 
catalytic efficiency (k cat /K m ) of the MERS-CoV M pro M298R mutant is higher than that of 
SARS-CoV M pr ° (Table 2). 

Conclusion 

The crystal structure of authentic N-terminus MERS-CoV M pro was determined and this was 
found to involve a dimeric form with less intermolecular polar interactions. Biochemical and 
AUC studies indicated that MERS-CoV M pro shows almost the same proteolytic activity as 
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Fig 5. ITC titration of peptide substrate binding to MERS-CoV M pr ° wild-type (A), E169A (B) and M298R 
(C) mutants. Raw data are shown in the top panel of each panel and represent the power input into the 
sample ampoule over time. Integrated data are shown in the bottom panels in terms of the total energy 
required for equilibration as a function of the molar ratio of substrate to M pro . The solid circles show the 
observed values and the lines represent the fitted results by ligand binding analysis using the Digitam 
program (TA instruments). The best-fit parameters for stoichiometry (N) and K d are shown in the bottom 
panels. 


doi:10.1371/journal, pone.0144865.g005 


SARS-CoV M pr °; although it is a monomer in aqueous buffer and displays substrate-induced 
dimerization (Fig 6). A conserved residue, Glul69, plays an essential role in the substrate- 
induced dimerization of both MERS-CoV M pro and SARS-CoV M pro . Moreover, mutation of a 
residue in the dimerization interface, M298R, was found to result in a more stable dimer form 
in aqueous buffer that had higher enzyme activity; while other two mutations, V4R and T126S, 
showed the reverse effect. Critical assessment of the residues important to dimerization of and 
catalysis by MERS-CoV M pro provides valuable insights into the mechanism that controls the 
monomer-dimer switch of important and valuable enzyme. 



Subunit 
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MERS-CoV M pro 


Subunit 

Association 


s M298R mutant, 

^ >\SARS-CoV M pro like 

Catalysis 



ES complex Catalytically Competent 


Fig 6. Schematic model of the catalytic process of MERS-CoV M pro . Substrate binding results in an 
active protomer (square) and then dimer formation. After a catalytic cycle, the dimer will dissociate to form the 
inactive protomers (circle). The M298R mutant, like SARS-CoV M pro , is able to maintain a dimer form without 
dissociation before the next round's substrate binding and catalysis. 

doi:10.1371/journal, pone. 0144865.g006 
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Supporting Information 

51 Fig. Sequence alignment of betacoronaviral M pro . Modified from an output from ESPript 
[36]. The green ovals indicate the catalytic dyad, while the blue ovals indicate the residues mak¬ 
ing intermolecular polar contact in the dimer interface of SARS-CoV Mpro. Magenta oval 
indicate the residue Glu playing dual role for the dimer interface and substrate binding site. 
Accession numbers are as follows: MERS-CoV, NC_019843.2; Bat-CoV_HKU4, ABN_010865.1; 
HCoV_HKUl, NC_006577.2; MHV_A59, NP_068668.2; SARS-CoV, NP_828863.1. 

(PDF) 

52 Fig. Expression and purification of recombinant MERS-CoV M pro . Protein identification 
by SDS-PAGE. M: molecular marker. Lane 1-5: cytoplasmic fraction, flow-through, elute from 
the nickel affinity column, flow-through after 4-h’s PL pro treatment and protein fraction from 
S-300 gel-filtration column. 

(PDF) 

53 Fig. Comparison with the structure of MERS-CoV M pro with that of SARS-CoV M pro 

(A) and those of ligand-bound complex, dimeric C148A mutant and bat-CoV HKU4 M pro 

(B) . (A) An overlay of the current structure of MERS-CoV M pro (cyan and orange) with that of 
SARS-CoV (grey; PDB entry luk3). The red arrows show the orientation change affecting the 
two domain Ills. Spheres show the two catalytic dyads. (B) Overlay of the current structure 
with ligand-bound complex (red; PDB entry 4YLU), dimeric C148A mutant (magenta; PDB 
entry 4WME) and bat-CoV HKU4 M pro (yellow; PDB entry 2YNA). 

(PDF) 

54 Fig. Activity assay of SARS-CoV M pro . The plot of rate constant (k 0 |,J versus the concen¬ 
tration of TQ6-pNA are indicated. The line represented the best-fit results according to the 
Michaelis-Menten equation (Eq 1). The protein concentration was 1.1 pM. The assays were 
performed in 10 mM phosphate (pH7.6) and repeated twice to ensure reproducibility and the 
error bars were shown. The kinetic parameters are shown in Table 2. 

(PDF) 

55 Fig. Continuous size distribution of MERS-CoV M pro at various protein concentrations. 

The protein concentration of 1.5 (solid line), 6 (dotted line) and 15 pM (dashed line) of 
MERS-CoV M pro were used and monitored the size distribution by AUC. The best-fit results 
suggest that the major species was a monomer with minor shift of sedimentation coefficient 
(1.74 to 1.97), while the calculated molar mass was from 35.8 to 36.6 kDa. The residual bitmaps 
of the raw data and the best-fit results are shown in the insets. 

(PDF) 

56 Fig. Effect of substrate concentration on the dimerization of MERS-CoV M pro . (A) Con¬ 
tinuous c(s) distribution of the enzyme at peptidyl substrate (TQ6-pNA) concentrations of 

0 pM (solid circles), 80 pM (open circles), 200 pM (solid triangles), 340 pM (open triangles), 
400 pM (solid squares) and 600 pM (open squares). The protein concentration was 0.25 mg/ 
ml. (B) Sedimentation coefficient shifts of the major species of MERS-CoV M pro at different 
TQ6-pNA concentrations. 

(PDF) 
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